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The magnetic phase diagram of the quasi one-dimensional spinchain system 
CsNiF3 below the Neel temperature is determined. For magnetic fields perpendicu- 
lar to the spin chains incommensurate phases are predicted. From linear spin-wave 
theory we obtain the instability line of the paramagnetic phase as a function of the 
strength and the direction of the field. The system undergoes a transition to a com- 
mensurate or an incommensurate phase depending on the direction of the magnetic 
field. In the commensurate phase the characterizing wave vector is locked to values 
describing a two-sublattice structure, whereas in the incommensurate phase the wave 
vector changes continuously between the corresponding two-sublattice wave vectors. 



PACS numbers: 64.70 R, 75.10 J, 75.30 D, 75.50 E 



Typeset Using REVTEX 



The quasi-one-dimensional magnetic compounds ABX3 (A alkaline metal, B transition 
metal, X halogen) have attracted much interest Jl|J^]. These materials have ferromagnetic or 
antiferromagnetic intrachain interactions and weak antiferromagnetic interchain interactions 
in the plane. The systems with integer spin (e.g. CsNiCla) are studied in context with the 
Haldane conjecture ||. CsFeCl 3 is discussed as an example of a system with a singlet- 
ground-state [[[]. 

CsNiFs is an example of a ferromagnet with planar anisotropy which has been studied 



extensively experimentally and theoretically |§-f|]. It remains in the center of focus fllOlJlT 
because it is a model system with reduced dimensionality where enhanced fluctuations have 
a pronounced effect on the ordering structure. Due to the planar anisotropy the spins will 
be oriented in the plane perpendicular to the chain axis. Neutron scattering revealed the 
dynamics of linear and nonlinear excitations ||. The one dimensional spin chain has no 
long range order but when applying a homogeneous field perpendicular to the chain axis 
spin waves can be measured and described within a linear spin wave theory || . However, the 
three-dimensional properties of CsNiF 3 have hardly been studied. Recently Baehr et al. [TI| 
measured the magnetic excitations in the three dimensional ordered state (T < = 2.7K). 
It could be shown that an isotropic antiferromagnetic exchange in the hexagonal plane and 
the dipole-dipole interaction are responsible for the three dimensional collinear long range 
order. The spins are oriented within the plane perpendicular to the chain axis along the 
crystal axes of the triangular lattice. 

In the following we study the three dimensional magnetic structure of CsNiF 3 for mag- 
netic fields oriented in the plane. We examine the stability of the paramagnetic phase as 
a function of the direction of the field. By means of linear spin wave theory we obtain 
the instability line at which the paramagnetic spin orientation gets unstable and changes 
to a canted structure. It turns out that the ground state depends sensitively on the field 
direction. For certain angular domains (A<£>j C pa 15.6°) there exist incommensurate phases 
separated by commensurate phases (of width Aip c pa 44.4°). Thus as a function of the direc- 
tion of the field we obtain alternating commensurate and incommensurate spin structures. 



This is the first time that incommensurate phases have been predicted in CSN1F3. 
The starting point of our investigation is the Spin Hamiltonian 

h = -2 j E s,s m + a £ (s; f -EE (W* + A w) s?sl - ^h £ s, . (i) 

i I a,/3 1,1' I 

Here J denotes the ferromagnetic nearest-neighbor intrachain interaction, A the single-ion 
anisotropy, J' w the interchain and i indicates positions on one and the same spin chain, 
whereas / indicates all spin positions. H is an external field perpendicular to the chain 
axis. The magnetic lattice structure of CSN1F3 is a simple hexagonal structure with lattice 
constants c = 2.6 A of the spin chains and a = 6.2 A of the triangular lattice in the plane. 
From neutron scattering |10| the coefficients in Eq. (P are deduced to J = 11.8 K, A = 3.3 



K and J' = —0.025 K. The value for J' is of the same order as the dipole energy which is 
necessary for stabilizing a collinear antiferromagnetic spin structure. Because of the large 
planar anisotropy the spins are forced to lie in the hexagonal plane. The classical ground 
state for vanishing field is given by three possible domains A-C (s. Fig. |I]) in which a 
collinear antiferomagnetic spin structure is realized, i.e. the rotation symmetry is broken 
due to the competition of the exchange and the dipole-dipole interaction |ITl[T2j] . Thus the 



system has an easy axis anisotropy. In the real system all three domains are simultanously 
present. 

Fourier transformation of the Hamiltonian (Eq. ([!])) yields 

# = -EE ( <V Q/3 - A5 az 5? z + J'J * + Af ) - g L fi B VNK S , (2) 

a,/3 q 

with the nearest-neighbor exchange energies (intrachain and interchain) 

J q = 2 J cos cq z (3) 
V o V ( ,0 aqx V3aq y \ 

J CL = 2J I cos aq x + 2 cos — - cos — — - I . (4) 

A^ denotes the Fourier transform of the long-range dipole interaction, which is calculated 
by means of the Ewald summation technique [ |I~3| , |I~4Tj . 

The qualitative effect of a homogeneous magnetic field is as follows. For weak magnetic 
fields transverse to the spin orientation the spins will reorient in order to gain energy from 



the Zeeman term. Thus when orienting the field along one domain direction the two others 
will change their ground state immediately. For a field longitudinal to the spins a spin-wave 
calculation reveals lT2| , [T5|j that the Neel state is at least metastable up to a finite critical field 
which depends solely on the dipole energy. Thus the situation for a virgin probe neglecting 
effects from domain wall energies and crystal defects is as follows: Without magnetic field 
the system might be built up of the three domains in equal fractions. Raising the magnetic 
field (parallel to the spins in domain A) does not change the spin orientation in domain A 
but leads to a slight reorientation in domains B and C. Above the the critical value the spins 
in domain A flip (first order phase transition) to an orientation identical to either domain B 
or domain C. For strong magnetic fields one finally enters the paramagnetic phase. When 
the magnetic field is decreased thereafter, the spins order again in the two domains but 
domain A is not formed any more because of the metastability of this domain. The system 
ends up in a state, where only domains B and C are present. 

In the following we study the excitations by means of linear spin wave theory. Via a 
Holstein-Primakoff-transformation we transform the Spin operators Sf in the Hamiltonian 
(Eq. to Bose operators a; and a\ JT5| , P~T| ] and diagonalize the quadratic form. 



In the paramagnetic phase all spins are aligned along the magnetic field H . Because the 
ground state of CsNiF3 for vanishing fields is not invariant with respect to a rotation around 
the spin-chain axis (recall that there are three domains A-C), the direction of the field plays 
a crucial role. The field direction is parametrized by </?, which is the angle between the field 
and the x-axis (inset in Fig. (4)). To quadratic order in the Bose operators we obtain the 
Hamiltonian 

H = Ep M + Y, A d + 9 ( fi q fl q G -Q + 5 q a q°-q) ( 5 ) 

q 1 

with the coefficients 

A q = S{2{ J - J q ) + A + 2J - 2 jy + g L ti B Ho 

+S (cos 2 <p(2Al x - A™) + sin 2 (p(2A v y - A™) - A zz - sin 2<pA%) (6) 
fi q = S(-A + A zz - sin 2 <pA x Q x - cos 2 <pA™ + sin 2ipA%) 



and the classical ground state energy of the paramagnetic state 

Ef M = -NS 2 ( J + J + (cos 2 <pA™ + sin 2 V Af)) - g L fx B NSH , (7) 

which is independent of the direction of the field for spherical shaped systems ( Aq x = Aq V ) . 
Here we considered only wave vectors within the plane (q z = 0, A y * = 0) due to the strong 
planar anisotropy. The dispersion relation for the paramagnetic phase is calculated via a 
Bogoliubov transformation to 

E* = ^Al-\B^. (8) 

This equation is valid for high magnetic fields. When lowering the magnetic field the para- 
magnetic phase gets unstable and changes to a canted spin structure. This instability is 
signalled by a soft mode at the wave vector q(ip) characterizing the phase below the para- 
magnetic phase. The value of the critical field Hq and the wave vector q(</?) depend on the 
angle of the field (p and can be evaluated from Eq. (|8|) by setting the excitation energy to 
zero. This leads to the equation 

g L ^ B H^) = 2S ( j; - J' + sin 2 <p{A% - A% y ) + cos 2 ip{A™ - A») - sin 2^) . (9) 

Note that this expression is independent of the ferromagnetic exchange and the anisotropy 
energies as long as they are much larger than the antiferromagnetic exchange and the dipole 
energies. Regarded as a function of q the maximum of this expression gives the angular 
dependent critical value H^ip) and the wave vector q(y?). Assuming that the phase joining 
the paramagnetic is a conventional spin- flop phase (a = (3 in inset of Fig. (4)) or a general 
two-sublattice structure, the so-called intermediate phase (a ^ (3 in inset of Fig. (4)), we 
expect the maximum value for the critical field at wave vectors describing the antiferromag- 
netic domains, i.e. qi = ^(0,1,0), q 2 = f (1,1/VM) or q 3 = f (1,-1/^3,0) (s. Fig. 
2). However, the detailed analysis shows that the paramagnetic phase gets unstable at an 
incommensurate wave vectorq(^) for certain field directions. Before considering the general 
case we study the special angles = 0° and <p = 90°. 
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(i) ip = 90°: For magnetic fields parallel to the y-axis the critical value is given by 

g L H B H? = 2S( j; - J' - Af + A?) . (10) 

This expression is found to have its maximum value of at qi, the antiferromagnetic wave 
vector of the collinear phase of domain A. The actual value for the critical field is evaluated 
with the parameters given above to (spherical shape) 

# cs/ = 340mT at q(90°) = q* . (11) 

At Hq V the system undergoes a transition into a commensurate phase, precisely to the two 
sublattice phase described by qi. Owing to the hexagonal symmetry, the critical field for 
ip = 30° is the same as for ip = 90° but at the antiferromagnetic wave vector q 2 characterizing 
domain C; for ip = —30° the structure is given by q3 corresponding to domain B (s. Fig. 
(5))- 

(ii) ip = 0°: For this field direction Eq. (§) reduces to: 

g Lf i B H™ = 2S(j; - J' - Aq x + A") . (12) 

The maximum value is not achieved for any of the two-sublattice wave vectors but for a wave 
vector with only an q x .-component depending on the relative strength of the antiferromag- 
netic exchange and the dipole energy. This follows from the fact that the dipole component 



A y y has a linear wave vector dependence at qo rather than a quadratic as found for J q ||18 
Evaluation of the critical value for CsNiF 3 leads to (spherical shape) 

= 290mT at q(0) = -(1.023, 0, 0) . (13) 

Thus the system undergoes a transition to an incommensurate phase. The incommensurate 
wave vector happens to be near the wave vector q 4 = ^(1,0,0), characterizing a four- 
sublattice structure, which describes an antiferromagnetic modulation along the x-axis. 

(iii) arbitrary ip =: Finally we turn to arbitrary angles, for which the situation turns 
out to be nontrivial. The complete dependence on the field direction q(<£>) for CsNiF 3 is 
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plotted in Fig. (|3]). The wave vector components q x and q y are plotted as a function of 
the angle ip for the region in question. A critical angle ip c ~ 7.8° is obtained above which 
the wave vector is locked to q2. Varying the angle (p between 0° and 7.8° the wave vector 
where the instability sets in changes continuously from q(0) to q(<£> c ) = q2- In Fig. (|J) 
the critical value for the magnitude of the magnetic field is plotted as a function of the 
angle for the same angular domain. Note that the critical field is continuous even at the 
critical angle ip c . It has only a small kink at this point. The dashed curve results from 
Eq. (g) under the assumption that the instability point occurs at q2 for the whole angular 
segment. Herefrom we can see that for ip < (p c the incommensurate structure is favored. Due 
to the inversion symmetry of the lattice together with q(ip) there is a second modulation 
wave vector — q(y?). The result of this investigation for all angles ip of the magnetic field is 
summarized in Fig. fl5|): The instability line of the paramagnetic phase is shown. There are 
angular regions where the paramagnetic structure undergoes a transition to a two-sublattice 
structure (annotated by the corresponding wave vector) seperated by regions drawn in thick 
where an incommensurate structure is formed. 

Examination of the commensurate phase leads to the following results: The commensu- 
rate phases are described near the transition to the paramagnetic phase by a two-sublattice 
structure. Calculation of the classical ground state energy |12[ shows that a conventional 
spin-flop phase which is parametrized by a single angle for both sublattice spins is stable 
only for the field directions of (p = 30° + n60° and integer values of n. For all other field 
directions the commensurate phase is an intermediate phase with two independent angles. 

In summary we have studied the ground state for CsNiF 3 with a homogeneous magnetic 
field oriented in the hexagonal plane. By stability investigations of the paramagnetic phase 
we obtained a non-circular instability line for fields in the hexagonal plane. The magnitude of 
the critical field and the type of phase joining the paramagnetic phase depends crucially on its 
direction. For certain angular domains the system changes to incommensurate structures, 
which are seperated by commensurate (two-sublattice) ones (Fig. ([])). This is the first 
time a staircase like behavior is found in a magnetic system as a function of the direction 
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of the field. The number of steps is 6, which corresponds to twice the number of two- 
sublattice wave vectors. This staircase is quite different from a devil's staircase with an 
infinite number of lock-in steps . Our results for CSMF3 are in contrast to Yamazaki 



et al. ||, who predicted conventional spin-flop phases for all field directions. This may 
result from their semiclassical model which does not consider the full nature of the dipole- 
dipole interaction. According to our theory the magnetic phase diagram shows a much 
richer structure including intermediate spin configurations and incommensurate phases. The 
observation of our predicted new phases and the wave vector dependence is left to future 
experiments, e.g. neutron scattering. 
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FIGURES 

FIG. 1. The ground state for CSN1F3 in the hexagonal plane is one of the three shown con- 
figurations (A-C), called domains. In domain A the two primitive vectors are represented. The 
antiferromagnetic modulation can be described by qi, q2 and q3 for domain A, B and C respec- 
tively. 

FIG. 2. The Brillouin zones of the hexagonal plane. The hexagon is the crystallographic and 
the rectangle (dashed) the magnetic one. 

FIG. 3. Wave vector (q x ,qy) at which the paramagnetic phase gets unstable versus the di- 
rection of the external magnetic field (solid line: q y , dashed line: q x ). Parameters for CSN1F3 are 
used. For ip > 7.8° the instability appears for q2. 

FIG. 4. Critical field below which the paramagnetic phase gets unstable for field direction 
between 0° and 30°. For smaller (larger) angles the paramagnetic phase changes to an incommen- 
surate (commensurate) phase. The dashed curve indicates the critical field when assuming that 
the soft mode occurs at q2- The inset shows the coordinate system for a general two-sublattice 
spin orientation. 

FIG. 5. Angular dependence of the instability of the paramagnetic phase for CSN1F3. 
The thick segments on the instability curve correspond to directions of the magnetic field 
Ho = Ho (cos <p, sin ip, 0) for which the paramagnetic phase changes to an incommensurate phase. 
These alternate with segments where a transition to a commensurate structure appears. The 
corresponding wave vector is given. 
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